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Abstract We show examples of excitation of coronal waves by flare-related 
abrupt eruptions of magnetic rope structures. The waves presumably rapidly 
steepened into shocks and freely propagated afterwards like decelerating blast 
waves that showed up as Moreton waves and EUV waves. We propose a simple 
quantitative description for such shock waves to reconcile their observed propa- 
gation with drift rates of metric type II bursts and kinematics of leading edges 
of coronal mass ejections (CMEs). Taking account of different plasma density 
falloffs for propagation of a wave up and along the solar surface, we demonstrate 
a close correspondence between drift rates of type II bursts and speeds of EUV 
waves, Moreton waves, and CMEs observed in a few known events. 
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1. Introduction 

Some solar eruptions are accompanied by large-scale wave-like disturbances vis- 
ible in various spectral ranges. Moreton waves (Moreton, 1960) observed in the 
Ha line have been initially proposed by Uchida (1968) to be a chromospheric trail 
of a coronal fast-mode MHD wave. Observations of the low corona in extreme 
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ultraviolet (EUV) with EIT (Dclaboudiniere et ai, 1995) on SOHO revealed 
large-scale wave-like disturbances visible as fronts of enhanced (but still low) 
brightness, either quasi-stationary or propagating over large distances up to the 
whole disk along the solar surface or expanding above the limb. These transients 
called "EIT waves" (or "EUV waves") are registered by a number of EUV 
telescopes — EIT, TRACE, STEREO/EUVI, SDO/AIA, etc. (see Warmuth, 
2007; Wills-Davey and Attrill, 2009; Gallagher and Long, 2010 for a review). 

Analysis of observations and interpretation of such phenomena meet prob- 
lems. Glaring flare emission hampers detection of faint EUV waves. Many data 
are limited by 12-min imaging rate of EIT. Faster TRACE observations (Handy 
et ai, 1999) have a small field of view. It is often difficult to reliably identify and 
trace a moving feature of interest. To overcome these difficulties, special methods 
are employed, but they might contribute artifacts. Multi-instrument analyses 
sometimes encounter timing problems. Thus, some results used in interpretation 
and modeling might not be completely adequate to the observed phenomena. 

The nature of EUV waves has been debated starting from their discovery 
(Thompson et al, 1998). The most popular interpretations of a near-surface 
EUV wave are i) traces of an MHD fast-mode wave in the lower corona {e.g., 
Thompson et ai, 1999; Warmuth et al. 2001, 2004b; Khan and Aurass, 2002; 
Hudson and Warmuth, 2004) and ii) plasma compression in bases of coronal 
loops in their successive stretching by an expanding CME {e.g., Delannee, 2000). 
A numerical 2D MHD simulation of a magnetic flux rope eruption (Chen, Fang, 
and Shibata, 2005) revealed both the {i) and {ii) disturbances. Schmidt and 
Ofman (2010) presented the first 3D MHD modeling of an "EIT wave" as a 
disturbance produced by an eruption-driven shock wave. A fast-mode wave de- 
tected in the simulation corresponded to an EUV wave observed in the modeled 
event including reflection from a coronal hole in support of the wave hypothesis. 

None of existing models describes all properties of EUV waves. Most likely, 
this is because the multitude of transients observed as EUV waves actually 
correspond to different phenomena. This conjecture is supported by a variety 
of morphologic and dynamic characteristics of observed EUV waves. For exam- 
ple, 1) their velocities estimated from observations of some events exceeded the 
coronal fast-mode speed, whereas they were lower in other events (Wills-Davey, 
DeForest, and Stenflo, 2007; Warmuth, 2010); 2) the wave front can be either 
diffuse or sharp; 3) kinematics of an EUV wave can be incompatible with the 
fast- mode MHD wave model (Zhukov, Rodriguez, and de Patoul, 2009). On the 
other hand, such properties of EUV waves as deceleration, decay, and broad- 
ening the disturbance (Warmuth et al. 2001, 2004a, 2004b; Long et ai, 2008; 
Veronig et al, 2010), bypassing regions of an increased Alfven velocity — coronal 
holes and active regions (Thompson et ai, 1999), possible reflections (Veronig, 
Temmer, and Vrsnak, 2008; Gopalswamy et ai, 2009) appear to correspond to 
the hypothesis of a coronal MHD wave. Considerations of wave-like transients 
sometimes observed in EUV to expand above the limb also suggest that different 
phenomena might be involved (see, e.g., Zhukov and Auchcre, 2004; Grechnev et 
al, 2006b, 2008; Pohjolainen, Hori, and Sakurai, 2008; Meshalkina et ai, 2009; 
Chertok, Grechnev, and Uralov, 2009; Cohen et ai, 2009). In paper HI (Grechnev 
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et al, 2011) we consider an event with a two-component "EIT wave"; a propa- 
gating component matched properties of a coronal shock wave, while a stationary 
component was presumably associated with a stretched CME structure. 

If some EUV waves are really due to coronal shock waves, then their corre- 
spondence with signatures of shocks in higher corona is expected. Type 11 radio 
bursts are considered as manifestations of shock waves propagating upwards in 
the corona {e.g., Vrsnak and Oliver, 2008). Klassen et al. (2000) concluded that 
almost all metric type II bursts were accompanied by EUV waves but stated 
no correlation between their speeds; on average, the speeds of EUV waves were 
three times less than estimates from drift rates of type II bursts. Biesecker et al. 
(2002) found that many EUV waves were not associated with type 11 bursts. 

Sheeley, Hakala, and Wang (2000) proposed that kinks or deflections of coro- 
nal rays at the flanks and rear ends of fast CMEs could be signatures of shock 
waves. An important common property of these flank/rear kinks was their decel- 
eration, while the authors did not reveal deceleration of CME leading edgcs^. A 
usually considered scenario for the formation of shock waves, which might show 
up in LASCO images, seems to be incompatible with wave signatures on the solar 
surface. In this scenario, when the CME speed exceeds the local fast-mode speed 
at heliocentric distances ^ 1.5i?0, a bow shock forms continuously pressed by a 
fast CME. However, a bow shock followed by a Mach cone can only be formed 
by a supersonic body of a fixed size, whereas CMEs expand omnidirectionally 
(with respect to their center). Hence, the conic bow shock geometry appears 
to be unlikely for wide CMEs. Since neither mechanisms nor heights of the 
shock formation have been established, possible association between surface EUV 
waves and CME components cannot be excluded. Indeed, Veronig et al. (2010) 
found the upper part of an expanding EUV dome to coincide with a white-light 
CME, while the lower skirt of the dome was a surface EUV wave. 

For these reasons it is difficult to expect that all observed properties of EUV 
waves could be explained by a single mechanism. Based on this assumption, 
our three companion papers are focused on those EUV waves, which are most 
likely associated with coronal shock waves. We address a few events, all of which 
were previously studied. Even for a rather uniform subset of phenomena prob- 
ably corresponding to similar processes, seemingly contradictions between some 
observational facts might occur. We endeavor to reconcile such discrepancies. 

One of challenging issues is reconciliation of EUV wave velocities and drift 
rates of type 11 bursts. This is a subject of the present paper I, the first one 
of the three companion papers. Assuming the shock wave nature of EUV waves 
under consideration, we try to settle disagreement between different studies. In 
this paper we use a simplest approximation of a self-similar shock wave that is 
convenient in comparisons with observations, which often show self-similarity of 
the wave front expansion. However, this approximation corresponds to the strong 
shock wave limit and cannot apply to all stages of events. A weaker shock seems 
to be more realistic, but its calculations are more complex. Propagation of a weak 



^The SOHO LASCO CME Catalog at http://cdaw.gsfc.nasa.gov/CME_list/ (afterwards 'the 
CME Catalog'; Yashiro et al., 2004) shows pronounced decelerations for three of the four 
events considered by the authors. 
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shock along the solar surface is analytically modeled in Paper II (Afanasyev and 
Uralov, 2011). Paper III analyzes propagation of a probable shock wave in the 
17 January 2010 event using both strong and weak shock approximations. 



2. Methodical Issues 

Excitation of shock waves in eruptive events seems to be undoubted, but their 
sources have not been established. Three possible exciters of shock waves are 
considered: i) a pressure pulse produced by a flare, ii) a super-Alfvcnic piston, 
and in) an impulsive piston. In case (ii), a bow shock continuously driven by a 
supersonic CME is expected to appear with kinematics determined by the driver. 
In case (i) or [iii), the shock wave expelled by an impulsive driver propagates af- 
terwards freely like a decelerating blast wave (c/. Pomoell, Vainio, and Kissmann, 
2008). Observations suggest that shock waves excited by impulsive drivers and 
freely propagating in the low corona do exist. This is indicated by deceleration 
of Moreton/EUV waves {e.g., Warmuth et al. 2001, 2004a), differences between 
propagation directions of the wave and a possible driver (Hudson et al., 2003), 
very early appearance of type II bursts. We accept this possibility as a working 
hypothesis and use a description of propagation of a blast shock wave. It is 
possible to calculate it analytically for two limits. One limit is a strong self- 
similar wave, whose length along a propagation direction is comparable with 
the curvature radius of the wave front. The opposite limit is a weak shock wave, 
whose length is much less than the curvature radius of the front and a typical size 
of inhomogeneities in the medium. Grechnev et al. (2008) found that a formal 
usage of expressions for propagation of a strong self-similar shock wave excited 
by a point-like explosion in a gas allowed to fit the speeds and positions of a 
Moreton wave as well as an EUV wave at the initial stage of the motion. 

2.1. Self-Similar Shock Wave Approximation 

Grechnev et al. (2008) used a simple model to describe propagation of such a 
blast-like wave in plasma with a radial power-law (PL) density falloff S from an 
eruption center, n — nQ{x/ho)~^ with x being the distance and no the density at 
a distance of Iiq. We use Iiq sa 100 Mm, close to the scale height. {Our notations 
are different from those used in the papers listed above). Self-similarity of a shock 
wave is ensured by a presumable large pressure excess inside the volume confined 
by the shock front over a non-disturbed medium. Propagation of a shock wave 
in the self-similar approximation is determined by plasma density distribution 
being almost insensitive to the magnetic fields. Such a wave decelerates when 
S < 3 due to a growing mass of swept-up material. Propagation of a strong shock 
vs. time t in plasma with a PL density model is described by an expression 

X{t) CX i2/(5-5) (1) 

The approximation has a singularity at a; — > (here also the wave velocity 
V — > oo); however, wave signatures are not observed at small distances. The 
approximation becomes inaccurate at large distances, being not limited from 
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below by the fast-mode speed. A wave traveling along the solar surface weakens 
at large distances and propagates, in the first approximation, in a flat-layered 
atmosphere. Expression (1) describes strong spherical shock waves, which seem 
to be unrealistic in solar conditions, but its usage within some range of distances 
can be justified. 1) An enhanced plasma density above an active region falls off 
both vertically and horizontally. A power-law description of the falloff seems 
to be acceptable. 2) The self-similar solution of a strong wave satisfactorily 
describes damping of a gas-dynamic shock wave up to Mach numbers M sa 2, 
when the wave is neither strong nor weak. 3) Applicability of gas-dynamic self- 
similar solutions to MHD blast shock waves is not obvious, because account of 
the magnetic field seems to be necessary. We note the following. With M ^ 2, 
the gas pressure behind the shock front exceeds the magnetic pressure, even if 
/3 = Cg/V^ <C 1 in non-disturbed plasma before the shock front; here M is a ratio 
of the shock speed to the fast-mode speed before the front, Cs and Va are the 
sound and Alfven speeds. That is, the plasma flow behind the shock front has a 
gas-dynamic character. The role of magnetic fields is also not crucial for medium- 
intensity shocks (Af > 2), which also strongly heat plasma, thus significantly 
increasing its pressure. For example, with M « 2, the plasma pressure behind 
the front of a wave perpendicular to the magnetic field is equal to the magnetic 
pressure before the front, even if /3 ^ 1 there. For a switch-on shock wave 
running along the magnetic field this occurs with a Mach number M w 1.5. 

Expression (1) was obtained under an assumption that the 5 index was inde- 
pendent of the wave propagation direction We will formally use this expression 
also when 5 = 5{'d), if variations of 5 are small with the change of the direction, 
i.e., ■ndS <^ di!). Note that in a limit of a weak, short shock wave, its propagation 
is determined by a local value of d even if this condition is not satisfied. The 
above considerations lead to a heuristic conclusion about a possibility to use 
expression (1) for approximate estimates of kinematic characteristics of shock 
waves of intermediate intensity propagating in medium with S = 6(19). 

It is useful to compare the power-law coronal density model with other pop- 
ular models. The Newkirk (1961) model {ue = 4.2 x lO'' x 10^-^^^'', r is the 
heliocentric distance expressed in solar radii) describes the radial plasma density 
distribution in a coronal streamer. The Saito model (Saito, 1970) describes the 
density distribution above the quiet Sun depending on the latitude (j) 

ne(r,(/)) 3.09,, „ . • 1-58,, „ . 0.0251, . n 5 /^n 

Figure la presents the Saito model for = 0°,60°, and the Newkirk model. 
The PL model can be adjusted to any of these models by varying its parameters. 
The riQ, Hq parameters are redundant; we have split them to clarify their physical 
meaning. The x variable in the radial direction is x ~ (r — l)i?0. The PL model 
with (5 = 2, no = 2.8 x 10^ cm^'^ agrees within ±30% with the Newkirk model at 
r = (1.2 — 9)i?o, which is important in anaylses of type II bursts. The parameters 
of the PL model can be adjusted to the Saito model for various (p as well, which 
is important in analyses of CMEs. A single PL model with a direction-dependent 
S provides a convenient alternative to complex involvement of various density 
models and allows one to account for individual properties of active regions as 
well as highly disturbed conditions just before the onset of a wave. 
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Figure 1. a) Coronal density models of Newkirk (dotted) and Saito (dashed, for latitudes 
= 0° and = 60°). The solid lines represent the power-law model with different parameters 
to fit the Newkirk model at r < 9Rq and the Saito model within the LASCO/C2 & C3 fields 
of view shown with shading. The vertical shaded region indicates the typical range of type II 
bursts, b) A cartoon illustrating the relation between a type II burst and an EUV wave. A 
narrowband type II emission is generated by a shock front propagating in a distinct extended 
narrow structure like a coronal ray. A near-surface EUV wave runs slower, (c) Speed -time 
plots of the EUV wave and the type II emission site. 



2.2. Shock Waves and Type II Bursts 

Assuming S to depend on a propagation direction, we get an approximation for 
a shock of an intermediate intensity propagating in an anisotropic medium (Fig- 
ure lb). A quasi- isotropic shock wave propagating in homogeneous corona can 
only cause drifting continuum radio emission. A strong narrowband harmonic 
type II burst can appear if the shock front passes along a lengthy structure 
like a coronal streamer {e.g., Reiner et al., 2003; Mancuso and Raymond, 2004). 
The cumulation effect increases the density jump in vicinities of the streamer's 
current sheet and intensifies radio emission (Uchida, 1974). The situation re- 
sembles a flare process running along a coronal ray (Uralova and Uralov, 1994). 
The difference between the horizontal and vertical directions explains a relation 
between the speeds of an EUV wave and a type II burst (Figure Ic). The upwards 
speed is higher than the surface one, because the wave center rises. 

Even if a shock wave appears at a zero height, its front rapidly becomes convex 
and tilted towards the solar surface (Figure lb). This front shape has been 
actually observed by Hudson et al. (2003) and is consistent with observations 
discussed by Warmuth et al. (2004b). If a convex shock front encounters a remote 
coronal ray, then the intersection site bifurcates, and its parts move along the ray 
in opposite directions (c/. Mancuso and Abbo, 2004). The contact corresponds to 
an infinite drift rate followed by bidirectional drifts to resemble a direct bracket 
"(" in a dynamic spectrum. Note that dynamic spectra present a combination 
of emissions originating at different sites, so that the intensities are summed. 

In Section 3 we reconcile kinematics of "FIT waves" and drift rates of cor- 
responding type II bursts in terms of our approach based on the self-similar 
shock approximation (Section 2.1). We apply power- law curves with the same 
onset time to both spectral domains (hereafter 'shock-PL fit'). The density falloff 
index in a streamer determining the drift rate of a type II burst is expected to be 
S ^ 2. Real lateral density falloffs in a streamer should be steeper than along its 
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axis. Otherwise, streamers would not be visible in homogeneous corona. Thus, 
real shock fronts should be oblate at medium distances. The density falloff for an 
EUV wave escaping from an active region can be < 5 < 3. For possible wave 
signatures in CMEs, the density falloff index is expected to be close to the Saito 
model (2), i.e., 6 ^ 2.6 at moderate latitudes and steeper at higher latitudes. 

2.3. Fit of Presumable Wave Signatures 

We fit the drift rate of a type II burst manually. The onset time to of a wave 
can be approximately estimated from observations. We take parameters of the 
plasma density model uq = 5.5 x 10* cm~^, Hq = 100 Mm, and 6 according 
to the considerations in the preceding Section. One more input parameter is a 
reference frequency /obs of a band with a harmonic number Nf (usually 1 or 
2) actually observed in a dynamic spectrum at a time ti. The corresponding 
plasma number density is ni = [fohs{ti) /Nf / {0.9 x IC*)]^, and the height is 
xi = ho {no/niY^^ . Then the height-time plot of the shock tracer is calculated 
as x{t) = xi [{t — to)/{ti — io)]^^^^^''^ the corresponding density variation as 
n{t) = no [x{t)/ho]^^ , and the outline of both bands of the type II burst as 
/(t) = [1,2] X 0.9 X 10'*[n(t)]^/^. By adjusting S and to in sequential attempts, 
we endeavor to approach a best outline of the bands. Uncertainties provided by 
the routine typically do not exceed 1 min for to and 0.2 for S. The spectrum can 
be coordinated with measured heights by adjusting no, as usually done. 

To fit presumable traces of shocks in coronagraph images (in this paper we use 
measurements from the CME Catalog), we employ two ways. The first way is a 
manual fit similar to the routine outlined in the preceding paragraph. Input pa- 
rameters are starting estimates of S and to, the heliocentric distances of the wave 
origin and the wave front ri measured at a time ti. The initial approximation 
of the height-time plot is r{t) = (ri - ro) [(t - to)/{ti - io)]^^^^~*^ + ''o- Then 
sequential attempts arc made to approach a best fit. The second way employs 
a log -log height -time plot, which is a straight line for a power law. We use a 
second-order fit and itcratively vary 6 and to to minimize the second-order term. 
One should be aware of the fact that the major but unknown uncertainties can 
be due to identification of the wave front in coronagraph images. 



3. Observations 

3.1. Event 1: 13 July 2004 

This event (Figure 2) was associated with an eruptive M6.7 flare (00:09-00:23, all 
times hereafter are UT) in active region 10646 (N13 W46) and a CME observed 
with SOHO/LASCO (Brueckner et al, 1995) after 00:54. Two parts of the CME 
(Figure 2b, c) are listed in the CME Catalog as two CMEs measured at position 
angles (PA) of 294° and 213°. A type II burst was recorded in three observatories. 
The three estimates of the shock speed progressively decreased in time. 

Grechnev et al. (2008) revealed signatures of a probable blast wave as an Ha 
Moreton wave (gray in Figure 2a) and an EUV wave (white). Both disturbances 
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Figure 2. The 13 July 2004 event, a) The Morcton wave (gray) and EUV wave (white) fronts 
superposed on the EIT 195 A difference image. The black arcs 1, 2, and 3 trace the directions of 
measurements. The dotted circle denotes the limb, (b) The leading part of the CME observed 
in two LASCO/C2 difference images at 00:54 and 01:54 separated with the black circle, c) The 
non-disturbed corona in a combined image composed from non-subtracted EIT 195 A, Mark4, 
and C2 images. The dashed ovals in (b) and (c) outline the CME in the 00:54 image. The 
dash-dotted straight lines in (b) and (c) mark the position angles at which the measurements 
listed in the CME Catalog were carried out for the fastest features Fl and F2. The axes show 
the distances from the solar disk center in arc seconds (a) and in solar radii (b and c). 



were kinematically close to each other and to a kincmatical curve expected for 
a lower trail of a decelerating coronal blast wave. Its exciter was not discussed. 
The authors proposed that the decreasing estimates of the shock speed reflected 
deceleration of a single shock wave, but they did not consider the type II burst. 

Pohjolaincn, Hori, and Sakurai (2008) [afterwards PHS] analyzed the type II 
burst in this event, but could not reconcile the overall drift with propagation of 
a single shock wave. The authors proposed that two shock waves were excited, 
one by a flare blast, and the second one by an expanding loop, a part of a CME. 

We carried out an additional analysis of this event, measured kinematics of 
an eruptive system in order to find out a probable origin of a shock wave(s), and 
to reconcile its (their) propagation with the EUV/Moreton waves and the CME. 

3.1.1. Eruptive System 

Figure 3 shows the eruptive system: a leading bright feature 'bf, two filament 
segments '1' and '2', and several eruptive loops, of which one ('loop') was con- 
spicuous and thus will be discussed henceforth. Long exposure times (33-46 s) 
caused a blurring of fast features, e.g., a jetlike appearance of the bright feature. 

To measure expansion of a feature, we outline it with an oval arc that allows us 
to trace it, even if its leading edge is difficult to detect. The obtained distance - 
time plot is used as a starting estimate. Then we choose a regular function to 
match the distance - time plot and estimate its parameters. Using the analytic fit, 
we calculate expected distance - time points, compare them with observations, 
and improve the fit. All kinematical plots arc calculated by means of integration 
or differentiation of the analytic fit rather than measurements. Our ultimate 
criterion is to reproduce the motion of an analyzed feature. Observational limi- 
tations do not allow us to reveal a detailed time profile of acceleration, and we 
describe it instead with a smooth bell-like function. We use a Gaussian time 
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Figure 3. The eruptive system in the 13 July 2004 event (TRACE, 173 A). The cross marks 
the initial position of the filament. Oval arcs outline the eruptive loop (dashed), filaments 
(solid) , and bright feature (short arc) . The broken lines denote the expansion directions of the 
loop (dashed) and the filament (dotted). Axes show arc seconds from the solar disk center. 



profile (sec Wang, Zhang, and Shen, 2009). Then acceleration a is 

a= K -Wo)exp{-[(t-<o)/Tacc]V2}/(V2^Tacc) (3) 

Here race \/8 In 2 is a full width at half- maximum of the acceleration time profile, 
which is centered at the to time; vq and vi are velocities at the onset and 
end of the acceleration stage. In cases of a more complex kinematics, we use 
a combination of Gaussians and adjust their parameters manually. 

Uncertainties are mainly determined by difficulties to identify and trace a 
feature in question. For this purpose we use both non-subtracted and difference 
images processed in various ways and improve results in sequential attempts. 
The worst traceable feature in this event is the eruptive filament, which appears 
as a semitransparent dark feature, whose segments 1 and 2 are faintly visible in 
Figure 3. We therefore present the measurements of the filament in more detail. 

The results of plane-of-sky measurements are shown in Figure 4 for the loop 
and bright feature (left) and for the filament segments (right). The dashed curve 
in Figure 4d was calculated from the dashed acceleration profile in Figure 4f to 
match an upper envelope of the measured positions of segment 1. The dashed 
curve in Figure 4e is a corresponding velocity profile. The dotted curves in Fig- 
ure 4d- f correspond to a lower envelope of the positions of segment 2. The black 
solid curves correspond to an average height - time plot. The horizontal bars show 
the exposure durations. We consider the maximum acceleration of the middle 
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Figure 4. Kincmatical plots of the eruptions in the 13 July 2004 event: the loop and bright 
feature (left) and the filament segments 1 and 2 (right; the solid eurves correspond to averages 
between 1 and 2). Symbols mark the measured plane-of-sliy distances, and curves represent 
their fit. The gray curve in panel (e) displays the plasma pressure computed from GOES 
SXR fluxes. The gray curves in panels (d-f) show kinematical plots calculated by assuming 
the correspondence of the acceleration plot to the microwave one. The vertical dashed lines 
delimit the start time of the wave estimated by Grechnev et al. (2008). 



part of the filament to reaeh about 2 km s^^ (—60%, +120%), and its probable 
speed and positions to be close to the solid black curves, being within the dotted 
and dashed boundaries. The bright feature and eruptive loop are visible much 
better, their leading edges are well defined, and therefore uncertainties of their 
positions are significantly less. The estimated height -time profiles match the 
measurements. The main uncertainty here is due to temporal undersampling, 
which is only crucial for the estimate of acceleration of the loop. 

When the measurements began, the bright feature already rose that probably 
corresponded to the initiation phase, which started at about 00:07 according 
to soft X-ray (SXR) GOES data. A strong acceleration started at about 00:13, 
reached « 4 km s~^, and then changed to deceleration. The loop was static by 
00:14:30; after 00:15:10 its speed sharply changed to Moop max ~ 320 km s~^ 
and did not increase afterwards. The last measured point hints at a possible 
later deceleration. The transition from the initial zero speed to a final one oc- 
curred between two samples. Hence, the maximum acceleration of the loop could 
well exceed 10 km s~^. All components of the corresponding CME decelerated 
(Grechnev et al., 2008), so that the eruption resembled an explosion with a 
strong impulsive acceleration followed by continuous decreasing deceleration. 

Figure 4c also shows the plasma pressure computed from SXR GOES fluxes 
and a source size of 15 Mm found from RHESSI (Lin et al., 2002) images. The 
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pressure gradually rose while the bright feature suddenly started to decelerate. 
Thus, the flare pressure was unlikely a driver of either the eruption or the wave, 
whose estimated start time is delimited with vertical dashed lines. 

The filament started to rise nearly simultaneously with the bright feature. 
However, both the acceleration and speed of the bright feature were higher, and 
it surpassed the filament (see Figure 3). The nature of this feature is difficult 
to identify. In some images it resembles an arcade surrounding the filament, but 
initially it seems to be located below the filament. This feature might be also 
one more filament, which brightened due to heating. An additional possibility is 
suggested by a scenario proposed by Meshalkina et al. (2009): this feature might 
be a small-scale magnetic rope whose eruption destabilized the filament. 

The gray curves in Figure 4d - f show the kinematical plots calculated under 
the assumption that the acceleration plot corresponded to the 9.4 GHz light 
curve (Nobeyama Radio Polarimeters, Torii et al., 1979; Nakajima et al., 1985; 
ftp: / /solar. nro.nao.ac.jp/ pub/norp/xdr/). Microwave bursts are known to be close 
in shape to hard X-ray (HXR) ones, while accelerations of eruptions have been 
found to be close to HXR bursts [e.g., Temmer et al., 2008). The plots calculated 
from the microwave burst lag behind the actual plots of the filament by about 
two minutes indicating that, most likely, the fiare was caused by the eruption. 
This fact suggests that the eruptive filament accelerated almost independently 
of the flare rcconncction rate and HXR emission, at least, in this event. 

There are two options regarding a relation between the bright feature and the 
loop. One possibility is that w 1.5 — 2 minutes after the start of the acceleration of 
the bright feature, the loop suddenly and independently underwent much higher 
impulsive acceleration. Alternatively, the loop was expelled by a shock front that 
appears to be more probable. For the latter case, the strength of the shock can 
be estimated. The Mach number is M = V^h/Vfast: where T4h is a shock speed, 
and T/fast is a fast-mode speed. The shock speed at the onset of the loop 
motion can be roughly estimated from a PL fit (Grechncv et al., 2008) to be 
about 1000 km s~^, but it is rather uncertain because of insufficient temporal 
coverage by TRACE images and their long exposure times. The fast-mode speed 
can be estimated from an expression Vs\i ~ Vfast + KC/sh/2, where k is a coefficient 
determining the steepening rate of the wave front. This coefficient, 1/2 < k < 
3/2, depends on plasma beta and the propagation direction. We take the speed of 
the loop as the gas speed behind the shock front, t/gh ~ V^oop max ~ 320 km s~^ 
(see Figure 4b). With these quantities and k sa 3/2 for the wave propagation 
perpendicular to the magnetic field in low-beta plasma, the Mach number is 
M w 1.3. One might suppose that the steepening time was about the interval 
between the peaks of the solid and dotted curves in Figure 4c corresponding to 
the acceleration profiles of the bright feature and the loop, i.e., about 0.5 min. 

We extrapolated coronal magnetic fields from a SOHO/MDI magnetogram 
using a package of Rudenko and Grechnev (1999) based on a potential approxi- 
mation (Rudenko, 2001). The result showed that the eruptive loops were strongly 
tilted to the photosphere, and their tops were at a height of about 30 Mm. With 
a height of the pre-eruptive filament of > 10 Mm, the shock front hit the loop 
almost horizontally, and the estimated Mach number is related to the horizontal 
direction, while in the vertical direction the shock was probably stronger. 
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Figure 5. a) Measurements from the CME Catalog for the Fl and F2 CME components fitted 
with the shoclc-PL model, <5 ^ 2.6. The circle denotes the position of the off-limb EUV wave 
at 294°. b) Speed — time plots corresponding to each fit. Linear-fit speeds are also specified. 
Shading denotes the interval in which the type II burst was observed. 



Thus, the loop itself was most likely passive, and its motion was driven by 
the shock wave. All observed products of the eruption monotonically deceler- 
ated starting from TRACE observations and up to LASC0/C3 ones. The loop 
therefore is unlikely to have excited the second shock wave, as PHS hypothesized. 

3.1.2. Moreton/EUV Wave, Type 11 burst, and CME 

Grechnev et al. (2008) could not find out from kinematics if the leading edge of 
a coronal transient observed by LASCO (see Figure 2b) was a mass ejection or a 
trace of a wave. Comparison with a non-subtracted image of the corona before the 
CME in Figure 2c suggests its spiky leading fringe to be coronal rays deflected by 
a wave (pronounced for feature Fl at PA = 294°). Vourlidas et al. (2003) and 
Magdalenic et al. (2008) interpreted such phenomena as signatures of shocks. 
Assuming that the CME leading edge was due to the wave, we outlined the 
measurements from the CME Catalog with a shock-PL fit. The calculated curves 
with S ~ 2.6 corresponding to the Saito model agree with the measurements 
(Figure 5). The decreasing speeds computed from these height -time plots at 
distances corresponding to the LASC0/C2 and C3 fields of view agree with the 
linear-fit speeds estimated in the CME Catalog, while within the interval when 
the type II burst was observed the speeds are higher by a factor of 2 — 2.5. 

Figure 6c shows a type II burst recorded by the HiRAS radio spectrograph 
and our shock-PL fit of both the fundamental and second-harmonic emission 
(solid lines). We use the density falloff index S ~ 2.1, which is close to the 
Newkirk model expected for a streamer. Features QFl and QF2 are open to 
question; their presence hinted PHS at two different shock waves. The dotted 
lines approximately reproduce the outline of PHS following the logic suggested 
by their Figs. 4 and 6. They correspond to fixed velocities of the type II exciters. 
However, a flare blast wave proposed by the authors is expected to decelerate: 
just our outline corresponds to a freely propagating blast wave. 
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Figure 6. The 13 July 2004 event. Height— time (a) and velocity - time (b) plots of the 
Morcton and EUV waves, and (c) the fit of the type II burst (5 = 2.1). Qfl and QF2 are 
features open to question. A reverse drift of Qf2 is detectable in the inset which shows a 
portion of the dynamic spectrum recorded in Learmonth with a higher spectral resolution. 



The solid shock-PL fit outhnes the whole slowly drifting structure from the 
decimetric range up to the lowest frequency. A question remains about features 
QFl and QF2. The former feature with an uncertain harmonic structure docs 
not seem to favor the dotted outline relative to the solid one. The weaker QF2 
feature, which PHS considered as the onset of the second type II burst, indeed 
seems to have a harmonic structure. Its shape in the higher-resolution spectrum 
recorded at the Learmonth station (US Air Force RSTN) resembling "(" outlined 
with a black arc suggests an encounter of a shock wave with a dense structure 
(see Section 2.2). These facts support association of the type II burst with a 
single decelerating shock wave. The drifting continuum, which PHS found to 
start at 00:13 (confirmed by the acceleration profile in Figure 4f), might be due 
to emission from outside of a prc-shock region expanding towards a decreasing 
density or, alternatively, from inside of the expanding region with a progressively 
depleting density, as PHS proposed. Compression of the environment in the pre- 
shock interval from 00:13 to about 00:14:50 might have produced an excessive 
plasma density, which we describe with a radial power-law falloff. 

For comparison we show in Figure 6a the distance - time measurements of the 
Moreton wave (triangles) and the EUV wave (open circles) from Grechnev et al. 
(2008) along with shock-PL plots for the three directions 1, 2, and 3 denoted 
in Figure 2a, and the dotted plot of a spherical weak shock propagating in 
uniform plasma calculated by using expressions of Uralova and Uralov (1994). 
Figure 6b presents the speeds for the strong (solid) and weak (dotted) shock 
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approximations along direction 1. Both approximations are close to each other 
far from the eruption center (r > i?i w 200 Mm) being somewhat different 
at shorter distances. The shock wave propagating along the surface probably 
became weak at r > when it left the active region and entered quiet Sun's 
areas where the coronal density and the fast-mode speed were nearly constant, 
i.e., S ^ 0. Closer to the eruption site (r < i?i), the plasma density presumably 
had a power-law falloff, and the shock was not weak. These conditions seem to 
favor the self-similar shock approximation, in which i) the shock propagation 
speed is proportional to its intensity and does not depend on the fast-mode 
speed, and ii) the wavelength is equal to the distance passed by the wave, i.e., 
its duration increases. The density falloff of (5 < 3 within an active region and at 
its periphery corresponds to deceleration and damping of such a shock. The limit 
of a strong shock is a convenient idealization to describe the formation stage of 
a single shock wave, which propagates far from its source. A real forming shock 
wave does not seem to be so strong that the decrease of the fast- mode speed could 
be neglected when the shock leaves the active region. This issue is beyond our 
scope. We only note that formation of the shock discontinuity in a disturbance 
produced by an impulsive piston presumably completes (and its intensity reaches 
maximum) in a region where the falloff of the fast- mode speed is steepest. 

Note that deceleration of the EUV wave sweeping over the quiet solar area was 
stronger towards the equator, as comparison of the three fronts in Figure 2a for 
the 1-3 directions shows. This is expected for a strong shock, whose deceleration 
is determined by the density distribution, which is maximum at the equator (see 
the Saito model). This is also expected for a weak shock, whose propagation is 
governed by the Alfven velocity decreasing towards the equator due to both the 
density distribution and the dipole magnetic field of the Sun. 

Our analysis of Event 1 has revealed a probable excitation of a single wave 
by an impulsively accelerated eruption and steepening into a shock within one 
minute. Then the wave freely propagated like a decelerating blast wave and 
probably formed the leading edge of the CME. Our results are consistent with 
the conclusions of PHS about the role of a rapidly expanding eruption, formation 
of the shock wave at a very low altitude, and their estimates of the shock speed. 

3.2. Event 2: 1 June 2002 

Meshalkina et al. (2009) revealed a possible coronal shock wave presumably 
excited by a collision of an eruptive magnetic rope with a magnetic obstacle in 
the 1 June 2002 event. An M1.5 flare (S19 E29) started at 03:50 and lasted only 
11 min. SOHO/EIT carried out the 'High cadence 195' program, and LASCO 
did not observe at that time. Figure 7 shows an off-limb EUV wave in this event. 

Figure 8a - f presents the eruption with extreme outlines (dotted and dashed) 
of its leading edge found from two sets of images processed in different ways. 
Figure 8g - i presents plane-of-sky measurements of kinematics using the same 
technique as for Event 1. The eruption accelerated up to ~ 7 km s~^ and 
then decelerated. The deceleration might be overestimated, because the eruption 
started to disintegrate and become transparent. Similarly to Event 1, accelera- 
tion occurred during the rise of the HXR burst recorded with RHESSI, while 
the plasma pressure gradually increased all the time. 



eit_waves_bursts_v2.tex; 19/01/2013; 22:28; p. 14 



"EIT waves" and Type II Radio Bursts 




b '^'i^^ ':::^yx^^^W c 




03:58:33 - 03:51 :38 04:05:33 - 03:58:33 04:1 2:34 - 04:05:33 



Figure 7. An EUV wave (outlined with ovals) on 1 June 2002 in EIT 195 A running-difference 
images. The dashed circles denote the solar limb. The slanted cross marks the flare site. 
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Figure 8. The eruption in TRACE 195 A images (a-f) and its kinematics (g — i). The dot- 
ted and dashed lines correspond to extreme outlines of the eruption. Panel (i) also shows a 
25 - 50 keV RHESSI time profile (gray) and the plasma pressure calculated from GOES data 
(dash-dotted). The vertical dotted lines mark the observation times of images a-f. 



Figure 9 shows kinematics of the EUV wave (a,b) and the type II burst (c) 
similarly to Figure 6. To reveal the harmonic structure of the burst, we use 
again the record made in Learmonth (the inset). The burst consisted of two 
pairs of emission bands with frequency ratios in pairs of 2.0 and « 1.5 between 
the pairs. The two pairs of bands resemble band-splitting usually interpreted 
as the plasma radiation from the regions upstream and downstream of a shock. 
However, the relative split is atypically large for the metric range (Vrsnak et at, 
2001). Alternatively, this situation suggests propagation of the shock front along 
two streamers located close to each other. We outline the burst structure with 
two harmonic pairs of power-laws 1 and 2 with 'f ' indicating the fundamental 
emission and 'h' the second harmonic. The difference in S (2.4 and 2.6) might be 
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Figure 9. The 1 June 2002 event. Height — time (a) and velocity — time (b) plane-of-sky plots 
of the EU V wave (black) and the fit of the type II burst source (blue) . Open circles in panel 
(a) represent the distances of the foremost wave fronts from the eruption center, and labels 
'a' - 'c' correspond to the images in Figure 7. Panel (c) shows the dynamic spectrum with a 
shock-PL fit. The vertical broken lines delimit the interval of the type II burst. 



due to differences of density falloffs in the streamers as well as different angles 
between the shock front and the axes of the streamers. The estimated wave onset 
time is 03:53:40, close to the acceleration peak time (Figure 8i). 

Figure 9a shows height - time plots corresponding to the drift of the type II 
burst along with the EUV wave expansion. The type II emission was observed 
when its source was presumably located at heights from 190 Mm (heliocentric 
distance of 1.27i?0) to 500-600 Mm [(1.7- 1.9)^©], i.e., lower than usually 
assumed. Moreover, the outline of the shock front matches the fast-drifting deci- 
metric continuum suggesting its relation to the shock wave, which presumably 
appeared still lower. The type II burst probably started when the shock front 
reached the streamer and ceased due to deceleration and damping of the wave. 

3.3. Event 3: 19 May 2007 

This event associated with a B9.5 flare at 12:48-13:19 (N07 W06) and a fast 
CME has been well studied due to efforts of several researches mostly from 
observations made with EUV Imager (EUVI) of SECCHI complex (Howard et 
al, 2008) on STEREO (Kaiser et al, 2008). Nevertheless, some questions remain. 

Long et al. (2008) measured kinematics of the EUV wave, found distinct 
deceleration, and stated that the low velocities of "EIT waves" could be due to 
their temporal undersampling. The observations were found to be consistent with 
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Figure 10. EUV wave observed with STEREO-A/EUVI at 171 A on 19 May 2007. A set of 
spatial profiles (b) computed from ratio images within the strip contoured in panel (a). The 
arrow points at the origin of measurements. The A, B, and C broken bars in both panels mark 
the reflection positions. The slanted broken lines in panel (b) outline the steepest slopes. 



an impulsively generated fast-mode magnctosonic wave or propagating MHD 
shock. However, they revealed an initial acceleration of the disturbance from a 
nearly zero speed. This mismatches an expected behavior of an MHD wave. 

Veronig, Temmer, and Vrsnak (2008) also found deceleration of this distur- 
bance indicative of a freely propagating MHD wave and revealed a wave reflection 
at a coronal hole. They assumed that the wave was initiated by the CME. because 
the associated flare was very weak and occurred too late to account for the wave 
initiation. They also revealed two eruptions following each other. 

Gopalswamy et al. (2009) measured propagation of reflected wave fronts and 
considered the reflections as an argument in favor of a wave nature of EUV 
transients (see also Schmidt and Ofman, 2010). However, Attrill (2010) proposed 
that the reflections resulted from a misinterpretation of the running difference 
data and suggested instead that two EUV wave fronts developed during the 
event. Indeed, running differences reliably show only the outer boundaries of 
expanding disturbances, while the inner picture reflects all changes occurring 
between two images subjected to subtraction (Chertok and Grechnev, 2005). 

To see whether or not the reflections actually occurred, we use a movie com- 
posed of non-subtracted 171 A images of STEREO-A/EUVI (euvi_ahead_171.gif 
in the electronic version of the paper) . A backward motion suggesting a reflection 
is visible after 13:11 northeast from a plage region denoted in Figure 10a. We 
repeated the measurements of Gopalswamy et al. (2009), but without any sub- 
traction. We only consider the flrst probable reflection in the direction exactly 
backwards, because it is difficult to reveal wave fronts in non-subtracted images 
of this complex event. Figure 10 corresponds to Figs. 1 and 2 from their paper. 
To enhance the sensitivity, we use in panel (b) spatial profiles computed as the 
sums over the width of each slice extracted from an image as panel (a) shows. 
Each image was normalized to a pre-event image at 12:16:30 (fixed-base ratios). 
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Figure 11. The first eruption on 19 May 2007 observed with STEREO- A/EUVI at 171 A. 
Top: non-subtracted images, bottom: fixed-base diff'erences. The arcs outfine the foremost edge 



of the eruptive loop system according to the fit shown with the soUd fine in Figure 10b. 



Figure 10b shows that the plage region inflected after arrival of the wave front 
and then returned back. A backwards motion from plage A is faintly visible 
after 13:14. Region B exhibited a weaker bend. It is not clear from which region 
of A and B was the wave reflected. The latter seems to be preferable, because 
continuations of both the direct and reflected slanted traces intersect farther from 
region A, while a prolonged standing of the wave at region A is doubtful. The 
slanted broken lines represent the speeds found by Gopalswamy et al. (2009) and 
agree with the slopes in Figure 10b. Thus, the results of the authors were correct, 
at least, for the first reflection. The wave reflected backwards was considerably 
slower than the incident wave. This fact supports the shock-wave nature of the 
disturbance. Indeed, if an incident shock wave propagating with a velocity Vine sh 
encounters a semitransparent 'wall' like a coronal hole, then the shock reflected 
backwards is slower: Vback sh ~ Vine sh — ^as, where V^as is a velocity of the gas 
trailing the incident shock front. Vgas can be up to the sound speed. 

What does the accelerating part prior to 12:50 display? Figure 11 shows EUVI 
171 A images with subtraction and without it. The outermost boundary of the 
expanding bright feature coincides with the edge of coronal loops visible in the 
earlier non-subtracted images. Then eruptive loops rapidly loose brightness due 
to expansion and become invisible in non-subtracted images. We conclude that 
the accelerating part measured by Long et al. (2008) was related to the expanding 
loops, while the decelerating part was related to the wave. It is difficult to 
distinguish an appearing wave, which brightens, from a piston (loops), which 
becomes transparent. For this reason the acceleration of the loops is uncertain 
within 120—270 m s~^; the outline in Figures 10 and 11 corresponds to 260 m s~^. 

Figure 12a shows the measurements of the wave presented by Veronig, Tem- 
mer, and Vrsnak (2008) and their shock-PL fit (thick blue curve). The wave start 
(12:50) corresponds to the early rise phase of the HXR burst (red, also from 
their paper). The power- law fit corresponds to the measurements of the authors 
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Figure 12. The 19 May 2007 event, a) Propagation of the EUV wave measured by Veronig, 
Temmer, and Vrsnak (2008, black) along with an HXR time profile (red), a derivative of the 
SXR flux (pink), and a shock-PL flt (blue), b) Type II burst with shock-PL-fitted four pairs 
of bands. Paired bands are shown with the same colors. See details in the text. 
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Figure 13. Coalescence of two shock waves (a) and height -time plots of CMEs observed on 
19 May 2007 (b) and 24 September 1997 (c). Symbols present data from the CME Catalog, 
lines show their shock-PL fit. The filled circle in panel (c) denotes the origin of the wave. 



better than their quadratic fit and the hnear one. The eruption in this event also 
accelerated before the appearance of manifestations of flare reconnection. 

The suggestion of Attrill (2010) about the second wave appears to be correct. 
The second eruption was probably triggered by the first one. Filament 2 activated 
at about 12:47 and erupted at 12:55- 12:57 according to TRACE 173 A images. 
The HXR time profile was complex, but two distinct episodes are detectable in 
the derivative of the SXR flux recorded with GOES (pink in Figure 12a). The 
onset times of the two waves were about 12:50 and 12:56. 

A shock wave trailing a preceding one must reach the leading front due to 
properties of shock waves; the two shock waves coalesce to produce a single shock 
front (and a weak backwards disturbance, which we are not interested in). Its 
speed is less than the sum of the initial fronts' speeds; however, the resulting 
shock is stronger and faster than cither of the initial ones. So the slope of its 
distance - time plot is steeper than the initial waves had, and its virtual onset 
time is later than for either of the initial waves, as Figure 13a outlines. 

Let us try to understand a complex dynamic spectrum in Figure 12b recorded 
in San Vito (USAF RSTN). A harmonic pair of weak bands If, Ih is sometimes 
detectable after 12:52. A strong type III burst at 12:55-12:58 probably associ- 
ated with the second eruption partially blocked the weaker bands If, Hi. Two 
stronger type II bands I'f, I'h appeared at 12:58. Their appearance looks like an 
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inverse-N-like shift of the initial bands to higher frequencies (thicker outhne in 
Figure 12b). The initial bands If, Ih are still detectable sometimes. The outline 
of both these pairs of bands has the same start time of the first wave, 12:50. The 
appearance of the second pair might be due to the entrance of a part of the shock 
front into a dense region located rather high above plage A (see Figure 10a,b). 
The surface EUV wave reached plage A slightly later, as expected for a convex 
front tilted towards the solar surface (see Section 2.2). 

The second shock front was probably manifest in weak bands 2f, 2h (pale, 
the start time is 12:56). The bands denoted (l+2)f and (l+2)h (orange outline) 
probably reveal the resulting shock with a virtual start time of 13:01. All of these 
bands overlap with others increasing the total emission at the intersections. We 
remind that various bands were most likely emitted from spatially different sites. 

We have not revealed manifestations of the two merging shock waves in im- 
ages or in the kinematical plot. This result could be expected due to different 
propagation conditions for a shock wave upwards and along the solar surface. 
The near-surface portion of a shock front decelerates stronger and experiences 
significant damping, thus becoming weak and approaching a linear disturbance. 
If the trailing shock wave succeeded to catch up or intersect with the leading 
one, the effect is expected to approach an interaction of linear disturbances, 
when two waves pass through each other experiencing interference, while the 
scheme in Figure 13a shows an essentially nonlinear effect. 

The CME was fast (958 km s~^) and decelerated. We assume that its leading 
edge was a trace of a wave and apply a shock-PL fit to the measurements in the 
CME Catalog at PA = 270° (Figure 13b). The onset corresponds to the virtual 
start time of the coupled shock wave, and the density falloff corresponds to the 
Saito model. A trailing poorly observed CME at PA = 310° and probably related 
to the same event was considerably slower (294 km s~^). Its deceleration might 
be due to the influence of the wave running ahead and problems of measurements. 

Our analysis does not pretend to be perfect, but it shows that even such a 
complex dynamic spectrum can be reconciled with EUV observations and CME 
expansion under assumption of the shock-wave nature of related disturbances. 
Oversimplified considerations of such a complex event can be misleading. In 
particular, Yang and Chen (2010) concluded that the EUV wave in this event 
ran slower in regions of stronger magnetic field that seemed to be a challenge for 
the wave hypothesis. However, 1) the authors considered the radial component 
of the magnetic field only, whereas the Alfven speed depends on its magnitude. 
2) By taking the range of the magnetic field strengths < 0.6 G, for which Yang 
and Chen obtained anticorrelation with the EUV wave speed, and a density of 
> 2 X 10* cm~^ from the Saito model, one obtains fi = {2/"f)C'^ /Vl Z 4, i.e., the 
wave must be almost insensitive to the magnetic field. 3) The fronts in their Fig. 4 
stretch west-southwest and become sharply pointed at 12:59, whereas Long et al. 
(2008), Vcronig, Temmer, and Vrsnak (2008), and Attrill (2010) showed the front 
to be blunt in this direction at that time. 4) Usage of the Huygens plotting to 
find trajectories of the wave front resulted in a strange picture of intermittently 
condensed and rarefied ray trajectories in their Fig. 4. Thus, the results of Yang 
and Chen (2010) do not offer problems for the shock-wave hypothesis. 



eit_waves_bursts_v2.tex; 19/01/2013; 22:28; p. 20 



"EIT waves" and Type II Radio Bursts 




-1000. 



-1000 



Figure 14. EUV wave (blue) and Moreton wave in the 24 September 1997 event, a) Non- 
subtracted EIT 195 A image with the first front (blue crosses). Blue dotted line outlines its 
foremost edge. Red fronts outline the Moreton wave. White broken lines trace the measure- 
ment great circles. Black crosses denote the origins of measurements. b,c) Difference ratio 
images with blue outlines of EUV wave fronts and red measurement great circles. The outlines 
correspond to Warmuth et al, 2004a (courtesy A. Warmuth). 



3.4. Event 4: 24 September 1997 

This event was associated with a short M5.9 flare (02:43-02:52, S31 E19). An 
Ha Moreton wave and EUV wave in this event were first analyzed by Thompson 
et al. (2000). Warmuth et al. (2004a, 2004b) found kinematical closeness of both 
wave fronts to each other and their deceleration. The first EUV wave front 
(Figure 14a) was sharp and bright suggesting that the main EUV-cmitting layer 
was low. A difference ratio image (Figure 14b) reveals weak wave manifestations 
south, southwest, and slightly west from the outline of Warmuth et al. (2004a). 
The third front in Figure 14c is close to their outline. Deceleration of the EUV 
wave was therefore even slightly stronger than the authors estimated. 

White and Thompson (2005) analyzed wave signatures in microwave images 
at 17 GHz but did not reveal any deceleration. They also found: i) the speed 
of the microwave disturbance was 830 km s~^ against « 500 km s~^ estimated 
for the Moreton wave; ii) the brightness temperature at 17 GHz was about five 
times higher than an estimate from EIT data, and the discrepancy could be 
reduced if the kinetic temperature at 17 GHz would be different (preferentially 
higher) from the characteristic temperature of the 195 A channel. The authors 
also concluded that the timing of images should be corrected by « 100 s for EIT 
and by ^ 180 s for Ha to reconcile all observations. These facts indicate that the 
layers emitting microwaves and EUV were not identical. The higher speed, lesser 
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Figure 15. a) Propagation of the EUV wave (squares) and Moreton wave (diamonds) mea- 
sured by Warmuth et al. (2004a). Filled gray circles show the of White and Thompson (2005) 
data, black open circles with error bars show them shifted by 40 Mm. Both data sets are 
outlined with shock-PL fit. b) The HiRAS dynamic spectrum outlined with shock-PL fit and 
a negatively drifting continuum (dotted), c) Presumable relation between the layers emitting 
microwaves and EUV/Ho. Thin vertical lines show cross sections of the wave fronts of the 
largest column emission measure contributing at 17 GHz. The star is the wave origin. 



deceleration, and higher brightness temperature (i.e., column emission measure) 
observed at 17 GHz with respect to EUV hint at a possibly higher location of 
the microwave-emitting layer. Figure 15 shows our suggestion in panel (c); panel 
(a) presents the distance - time plots from both papers. The corrected times of 
EIT images (-1-99 s) and Ha ones (-1-170 s) arc specified at data points. 

To reconcile kinematics of the microwave and EUV/Ha wave fronts, we shift 
the White and Thompson data by 40 Mm and fit both data sets with the same 
onset time of 02:46:50 but different power-law exponents. Figure 15c explains the 
idea: the lower part of the front propagating in high-density regions decelerated 
stronger (6 ~ 2.3 from Warmuth et al. measurements), while microwaves were 
dominated by long cross sections of the wave front (bars 1-4) running in lower- 
density regions, S ~ 2.8 (c/. Paper III). The large-height EUV wave's leading 
edge detectable close to the eruption center diminished at large distances, where 
the EUV wave was dominated by low structures. The wave presumably appeared 
at a considerable height (the star). With a difference between the origins of 
measurements (black crosses in Figure 14a) of sa 32 Mm, our shift of 40 Mm, 
and spherical wave fronts, we get a height of 117 Mm. It seems to be overesti- 
mated; an estimate of about « 75 Mm from the dynamic spectrum appears to 
be more plausible. Thus, the wave could be strongly anisotropic starting from 
its appearance, or the wave exciter was large, or both. 

The complex dynamic spectrum in Figure 15b has a questionable harmonic 
structure, which is beyond our scope; we are only interested in the drift of the 
burst envelope. We have outlined presumable fundamental and second-harmonic 
bands. A possible higher-frequency band and the highest-frequency envelope of 
the burst are formally outlined as the fourth and sixth harmonic (this does not 
mean reality of emissions at such high harmonics). The outline corresponds to 
the initial height of 75 Mm, the same wave start time of 02:46:50, and 6 « 2.1 
typical of streamers. A negatively drifting continuum (broken outline) at the 
initial stage indicates propagation of the shock front towards the chromosphere. 
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A poorly observed CME centered at 137° with a speed of 531 km s^^ was 
injected into a preceding CME. The CME Catalog estimates its acceleration to 
be positive with a remark about uncertainty. By adding the known origin of the 
wave (the filled circle in Figure 13c), we get a shock-PL fit of the measurements 
from the CME Catalog with an exponent (5 « 2.5 close to the Saito model. 

Our considerations confirm correctness of both Warmuth et al. (2004a) and 
White and Thompson (2005) results, reconcile them with each other and with 
the type II burst as well as the CME. Deceleration of the front portion detectable 
at 17 GHz was much less than its lowest part visible in EUV and Ha had. It was 
not possible to reveal deceleration from microwave observations, which allowed 
detection of the wave within an interval as short as 4.5 min. 



4. Discussion 

TRACE observations of abrupt eruptions in events 1 and 2 have revealed plane- 
of-sky accelerations of magnetic rope structures of 4— 7 km s^^, i.e., (15- 25)-fold 
gravity acceleration. Then the eruptions in both events decelerated. Coronal 
waves appeared in events 1-3 approximately at the peak of acceleration. The 
onsets of the waves in all four events corresponded to the rise phases of HXR 
or microwave bursts (in event 4 according to White and Thompson, 2005). The 
wave in event 1 steepened into a shock within one minute and reached a Mach 
number of about 1.3 in a horizontal direction, while the upwards shock could 
be stronger. In the next Section wc consider which of known shock formation 
scenarios (see, e.g., Vrsnak and Cliver, 2008) appear to match the observations. 

4.1. Comparison of Observations with Shock Formation Scenarios 

A scenario of the shock formation by a flare pressure pulse is based on an idea 
that the /3 ratio of the plasma pressure to the magnetic pressure sharply changes 
from /3 ^ 1 to /? > 1. Such a change of /3 is believed to be possible in a flare loop 
and considered as a cause of an omnidirectional disturbance. Dramatic changes 
of the volume of a loop or its abrupt motion are necessary to get a significant 
intensity of a wave excited in this way. Objections against this scenario do exist. 

1) As Grechnev et al. (2006c) showed, the effect of a high /3 (even /3 > 1) in 
a flare loop is not dramatic, only an increase of all its sizes of y/l + /3. 2) The 
idea that a situation of /? — > 1 can cause instability of a loop is not confirmed 
by observations. Ichimoto et al. (1993) and Grechnev et al. (2006a) showed that 
the < 1 condition was not satisfied in long-lived hot coronal loops. 3) The 
flare pressures in events 1 and 2 rose gradually, without any marks of the wave 
appearance. 4) RHESSI showed the SXR-emitting regions in events 1 and 2 
to be fixed when the pressure increased. Veronig, Temmer, and Vrsnak (2008) 
concluded that the wave ignition in event 3 by the flare was unlikely. Thus, 
the theoretical considerations supported by the observational facts make the 
ignition of waves by flares doubtful in agreement with the conclusion of Cliver 
et al. (2004) (note that the authors implied the shocks to be driven by outer 
surfaces of CMEs, whereas we consider the shocks to appear inside CMEs). 
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This conclusion is seemingly opposed by the results of Magdalenic et al. 
(2008, 2010) who presented five events, in which shocks appeared during the flare 
impulsive phase, whereas related CMEs were slow. However, the authors did not 
consider a rare type of CME kinematics with very strong initial acceleration 
followed by deceleration. An example is our event 1 (Section 3.1, Figure 4). The 
acceleration stage in this event lasted three minutes. The fastest feature accel- 
erated up to 470 km and continuously decelerated later on. From an overall 
height -time plot, including measurements from LASCO images, Grechnev et al. 
(2008) estimated the speed of this feature to approach 100 km s~^ far from the 
Sun. Nevertheless, the shock in this event was most likely excited just by this 
feature as an impulsive piston rather than by the flare. Thus, a slower motion 
of the main body of a CME well behind a shock front does not guarantee that 
the CME or its components were not implicated in excitation of the shock. 

The maximum plane-of-sky speeds of the eruptions in events 1-3 appear 
to be well below the Alfven speed expected at moderate heights (< 100 Mm) 
above active regions {e.g., Mann et al., 2003), where the waves appeared. It 
is possible that the angles between the velocity vectors and the line of sight 
significantly differed from 90°, so that the real velocities could be higher, but 
the corresponding corrections are insufficient to increase the velocities above the 
Alfven speed. The time profiles of the velocities estimated for the eruptions and 
waves were quite different. Thus, the bow shock scenario is also unlikely. 

4.2. Impulsive Piston Scenario 

In a simplest scheme, a piston moving with a speed U has a box-like acceleration 
profile with a value a during the acceleration phase and zero before and after- 
wards. An important condition is that plasma cannot flow around the piston (this 
occurs, e.g., in 3D expansion of an arcade). In a flat geometry and homogeneous 
medium the plasma flow ahead of the piston corresponds to a simple wave. 
The discontinuity appears at = Vfast/('^a) at a distance r^^ = l^ast^sh '^^^^ 
1/2 < «; < 3/2 that is similar to a solution of an analogous gas-dynamic problem 
(Landau and Lifshits, 1987). Then the speed jump in the discontinuity increases 
up to the piston's maximum speed t/max- The condition U > T4ast essential for 
bow shock formation is unnecessary in the impulsive piston scenario. 

Accelerations and their durations before the wave onsets were (Section 3): 
4 km s~^ and 90 s in event 1; 7 km s~^ and 70 s in event 2; 0.12 — 0.27 km s~^ 
and 800 s in event 3. The stronger acceleration, the faster a shock appeared. 

The shock waves in the four events were most likely excited by eruptive struc- 
tures as impulsive pistons, which one might call the appearing CMEs. Then the 
waves rapidly steepened into shocks, detached the pistons, and freely propagated 
afterwards like blast waves. The shock excitation mechanism implies a source 
height to be nonzero, but rather low, probably < 100 Mm, as suggested by all 
dynamic spectra and implied by event 4. Magdalenic et al. (2008, 2010) found 
the heights of sources of metric type II bursts to be between 70 and 280 Mm; 
however, the apparent heights of limb sources in the metric range could be 
noticeably reduced due to refraction (see, e.g., Zheleznyakov, 1970). The shock 
character of the waves is supported by the correspondence of their kinematics 
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to the expected propagation of shock waves as well as the drift rates of type II 
bursts and drifting continua. The shock- wave nature of EUV waves is supported 
with probable reflections and coupling of two shock waves in event 3. Expansion 
of leading edges of CMEs produced in three events corresponded to propagation 
of the lower skirts of the shock fronts observed as Moreton/EUV waves. Thus, 
the wave excitation by an impulsive piston appears to match all the considered 
observations, basically corresponding to a picture proposed by Uchida (1974). 

So far we did not relate a piston with a particular structure or its surface. 
Presumable pistons could be either an eruptive filament (EE) or a CME frontal 
structure (FS). Both expand as an entire ensemble in a completely formed CME; 
only its outer sheath can be a piston. This sheath is believed to be both the 
surface of contact discontinuity and the outer FS surface. Expansion of such 
an FS-piston determines propagation of an interplanetary piston-driven wave 
and the drag force affecting a CME. The situation is different during the early 
CME formation inside an active region, when EE moves faster than it would 
be adequate for a self-similar expansion of the whole CME. The EE acts here 
as an impulsive piston and excites inside a future CME a wave, which freely 
propagates outwards as a shock wave. In the four events we revealed just this 
excitation scenario of waves, which resembled blast ones. Propagating upwards, 
such a wave inevitably would pass through the FS and appear ahead of it. 

4.3. CME Components and Waves 

4-.3.1. Particularities of Expansion 

Expansion of a magnetoplasma CME's constituent is different from kinematics 
of a wave traced, e.g., with a leading edge of a plasma flow driven by a shock. The 
CME expansion is known to be about self-similar at moderate distances from 
the Sun. The self-similar approach does not apply to early stages of expansion, 
when the structure and shape of a CME have not yet been established. When an 
instability driving an eruption completes and drag of the solar wind is not yet 
significant, the self-similar CME kinematics can be obtained from considerations 
of forces affecting a CME (Low. 1982; Uralov, Grechnev, and Hudson, 2005). 

Expansion of magnetoplasma structures is governed by magnetic forces, plas- 
ma pressure, and gravity as long as the effect of the solar wind is small. With the 
polytropic index 7^4/3 all the forces integrated over the boundary and volume 
of a CME scale with distance r by the same factor of r^^. This leads to an ex- 
pression for the CME velocity V^^^^ = + {V^ ~ Y^){\ - Ro/r), Vcme = dr/dt 
(Grechnev et ai, 2008). Here Rq is the initial size of self-similar expansion, Vq 
the initial velocity at Rq, and Voo the asymptotic velocity in infinity. At large 
distances acceleration c>c — >■ and Vcme — > Voo - The expression for the CME 
velocity describes different types of kinematics. The situation Voo S> Vo ^ ap- 
pears to be typical. Event 1 showed a different behavior, Vq > Voo, resembling an 
explosion with an impulsive acceleration followed by deceleration. A special type 
Vq ~ Voo (an impulsive acceleration is required to reach Vq) might correspond 
to some jet-like ejections, where magnetic reconnection destroys a structure of 
an eruption (Meshalkina et ai, 2009; Filippov, Golub, and Koutchmy, 2009; Liu 
et ai, 2011). Event 2 might have belonged to this type. 
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Since a FS starts to expand practically from a static equilibrium in a typical 
situation Voo ^ Vq ~ 0, the FS-piston usually either accelerates or moves with 
a nearly constant speed at the self-similar stage. By contrast, shock waves in all 
considered events decelerated. Hence, a typical FS-piston is expected to even- 
tually approach the wave front. What does such a relation between the speeds 
of the shock front and piston mean? Do MHD equations allow a decelerating 
shock wave to run for a long time ahead of a non-decelerating FS-piston? We 
search answers in a theory developed by Low (1984) in solving a problem of 
self-similar expansion of a CME preceded by a strong shock wave. Though the 
problem was solved in a limit of very strong shock propagating in plasma with 
too steep density falloff r^^^/''', the solution correctly relates accelerations of the 
piston and piston-driven shock. Assuming a common linear profile of the plasma 
velocity in the whole region from the expansion center up to the shock front, it 
is possible to relate kinematics of the shock front and the contact discontinuity, 
i.e., the FS-piston without a complete solution of the problem. By fitting the 
motion of the piston with a function rpjgt = bt"^, we express the sign of the 
shock acceleration ash vs. to: ash oc {am — with a = (7 + l)/2. Thus, 

the conditions ash < and apist — <P r p^st / dt^ > are incompatible. Such an FS- 
piston and the shock front expand in different ways, which cannot be coordinated 
with each other, unlike a bow shock. An impulsively excited freely propagating 
shock wave must eventually change to a piston-driven mode. Presumably this 
typically occurs at large distances, probably beyond the LASC0/C3 field of view. 
The transformation of a blast shock wave into a piston-driven one marks switch- 
on the aerodynamic drag and termination of the self-similar expansion regime. 
The drag force becomes significant, which means establishment of a continuous 
energy transport from the FS-piston to the shock wave. By contrast, a blast-like 
wave excited by an EF and running ahead of FS, which does not yet act as 
a piston, facilitates expansion of a CME into the solar wind. The shock wave 
forwards a part of its energy to the FS-piston, and the drag force is absent. 

Most likely, real shock waves are neither purely blast waves nor purely piston 
ones. A shock front is sensitive to any events occurring behind it, e.g., changes 
of the FS-piston speed, because the fast-mode speed behind the shock front is 
higher than its phase speed. To produce one more shock wave, an FS-piston 
has to repeat the maneuver, which produced the first shock. This is improbable 
when an instability driving an eruption has completed and a CME left the Sun. 

4.3.2. Distinguishing between Shock Signatures and CME Components 

Sheeley, Hakala, and Wang (2000) and Vourlidas et al. (2003) considered distor- 
tions of coronal streamers as a morphological suggestion of presumable shocks. 
Indeed, moving wave-like kinks or deflections of coronal rays resemble an ex- 
pected effect of a propagating shock wave. Some fast "CMEs" in difference 
images might be actually combinations of coronal rays deflected by shock waves. 
However, Filippov and Srivastava (2010) demonstrated that deflection of coronal 
rays could be due to expansion of a CME in magnetized corona without a shock. 
Irrespective of a particular type of a possible shock, more reliable morphological 
suggestion might be a spike-like leading feature due to deviation of a coronal 
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ray by a wave. Magdalenic et al. (2008) showed such a situation in their Fig. 4, 
where deflected coronal rays were visible well ahead of trailing CME structures. 

We remind that the conic bow-shock shape is not expected, at least, for wide 
super- Alfvenic CMEs. The shock front must cling to its foremost edge and closest 
flanks, while far flanks and a rear part can be constituted by a freely propagating 
shock front, so that the shape of the whole front would resemble an egg. 

An attractive way to detect a shock front is to search for discontinuity in 
the density distribution shown by coronagraph images (Vourlidas et al., 2003). 
However, this way is model dependent. Besides modeling the coronal density, 
one has to distinguish between the shock discontinuity and the contact surface 
separating the CME and environment. The 3-dimensional contact surface should 
be also modeled. Thus, identifying a shock front in coronagraph images does 
not seem to be a simple task. Combinations of different indications seem to 
be be useful. A non-structured faint density enhancement forming the envelope 
of a transient could be among them (Vourlidas et ai, 2003). Fast decelerating 
halo CMEs with such edges might be shock candidates. A worthwhile shock 
indication might be such a halo edge crossing a distorted streamer. An important 
complement of morphological suggestions is kinematics of an expected shock 
wave and its correspondence to the drift of a type II burst. Paper III also address 
changes in shape of a shock front occurring in its propagation in the corona. 

4.4. Presumable Scenario 

Our observations and considerations suggest the following presumable scenario of 
a flarc-rclatcd eruptive event. An eruption occurs due to a rapid development of 
an instability in a magnetic structure. An abruptly accelerating eruption destroys 
a pre-existing magnetic configuration, thus causing a flare, and produces an 
MHD disturbance as an impulsive piston. The disturbance appears at a height 
of ^ 50 Mm during the rise phase of an HXR/microwave burst, leaves the piston, 
rapidly steepens into a shock, and then freely propagates like a blast wave. Its 
lower trail might be observed as a decelerating Moreton wave as well as an EUV 
wave, and the wave dome is sometimes observed to expand above the limb. 

The motion of the shock front shows up in radio spectra as a drifting contin- 
uum and, when the shock front reaches the current sheet of a coronal streamer, 
as a type II burst. For the fundamental emission this usually occurs at ^ 100 
MHz (r ~ 1.5i?0). Metric type II bursts are expected to cease due to damping of 
shock waves that typically occurs above 20 MHz (r < 3Rq). Revival of a shock is 
possible at a few Rq due to decreasing Alfven speed, and decametric/hectometric 
type II emission can appear. The complex piston-blast-piston transformations 
of shock waves traveling in the corona with significantly varying parameters and 
possible coupling of multiple shocks imply well-known disagreement between 
metric and interplanetary type II events {e.g., Cane and Erickson, 2005). 

Expanding shock fronts can form envelopes of CMEs. Measurements in the 
CME Catalog referring to a fastest feature might be related to shock waves 
for fast decelerating CMEs, especially halos. Since a shock wave decelerates, a 
trailing mass must eventually approach its front. The shock becomes a piston- 
driven one presumably at distances r > 20i?o. The aerodynamic drag becomes 
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important. This picture is consistent with results of several cited papers. The 
story of shock waves associated with flare-related CMEs appears to be more 
complex than often assumed, in fact combining different scenarios. 



5. Concluding remarks 

Our seemingly simplified approach has resulted in surprisingly fine reconciliation 
of EUV waves, Moreton waves, metric type II bursts, and leading edges of CMEs. 
The first consequence is that independent of the quality of our approximation, all 
these phenomena are really manifestations of a common agent, i.e., a traveling 
coronal shock wave excited by an eruption. Second, our approach indeed provides 
a promising instrument for analyses slow-drifting bursts and their comparison 
with other eruption-related phenomena. The power-law approximation turns out 
to work well beyond conditions, for which it was derived. Our results clarify 
relations between flares, traveling coronal shocks, CMEs, associated wave-like 
manifestations, type II bursts, and provide a common quantitative description 
for some of these phenomena. Our important by-product is an indication of the 
leading role of eruptions with respect to flares, i.e., that the acceleration of an 
eruption occurs almost independently of the flare reconnection rate. 

Our approach and analysis needs elaboration and continuation. A number of 
issues to be addressed remains. Low heights, at which type II emission sometimes 
appears, indicate that our consideration of its generation in the current sheet 
of a coronal streamer needs elaboration. Data sets similar to those analyzed in 
our paper should be compared with imaging observations in the metric range. 
The analysis of decimetric to metric drifting bursts should be extended to longer 
radio waves in conjunction with coronagraphic observations. Despite success of 
our self-similar shock approximation, a more realistic weak shock approximation 
should be considered. The last issue is a subject of our paper II. 
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